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Abstract 


Arc furnaces are probably one of the worst kinds of load on the power system. Arc 
behavior is often very random, giving rise to large amount of harmonics and a poor 
voltage profile at the Point of common coupling. Besides this, the system also suffers 
firom flickering of incandescent lamps at loads connected on the same bus. To overcome 
these problems an arc model has been structured so as to simulate an arc furnace 
installation. The system is compensated with S.V.C as well as controlled current sources. 
The results indicate that the current compensation scheme studied in this thesis gives 
considerably better results than the conventional compensation approach using S.V.C’s. 
The level of flicker has been reduced to a great extent (80-85%) , indicating the 
effectiveness of the scheme for fast compensation of flicker and reactive power. 



Chapter 1 


INTRODUCTION 


1.1 History of electric arc furnace steelmaking 

Formulations of steel and other metals have been known since Biblical times. When 
Tubal Cain scraped ashes from his fire and discovered that he had reduced Iron or 
smelted metal, it is probable that his first thought was to make a weapon of destruction. 

The electric furnace was invented by a Frenchman, Pichon for which he was granted a 
patent on March 16, 1853, for his claim of economically melting minerals and metals. 
The development of dynamo in 1867 made the electric furnace feasible [11]. 

The first practical application of using electric power for melting steel was by Sir 
William Siemens in 1879. He used a crucible with horizontal electrodes to form a single 
phase arc above the crucible. This principle is still in use despite its being inefficient and 





costly. Stassano, in Italy suggested an indirect arc. Rennerfelt proposed a design that 
consisted of, two inclined electrodes and a third vertical electrode so as to direct the arc 
towards the bath for increased efficiency [11]. 

The next step in the development involved direct arc furnaces. These designs were 
characterized by vertical electrodes and current travel from the electrodes to the bath . 
The main improvement was that the heat generated by the arc was held within the scrap 
during melting which tended to reduce erosion. 

The jSrst electric furnaces were usually small units of 1 to 15 tons. By 1930, production 
had increased to annual output of 686,111 net tons [11]. As new alloys were developed 
for tool, valve, and stainless steel applications the industry continued to expand. At the 
start of World War II, the technology had improved to the point where the electric 
furnace industry was ready to provide the new types and greater quantities of metal 
required for military purposes. 


1.2 Problems associated ¥i1th arc furnace operation 

During melting down of the cold charge within the furnace, the three phases act much 
like three single phases. Two electrodes can strike the charge and draw the current 
without any current flowing through the third electrode. Currents of varying magnitude 
flow in the whole system (including the supply system). These currents are limited only 
by the impedance of the arcs, which is essentially resistive, the impedance of the 
secondary circuit, which is partly resistive and mostly reactive, the impedance of the 
transformer and by the impedance of the supply system 

It is during the early part of the meltdown when these fluctuating currents cause the 
greatest disturbance in the supply system If not properly controlled they will cause 
objectionable flickering of lights that uses energy from the same power supply system 
[ 2 ]. 




As the heat in the fiirnace progresses, conditions improve. When the refining stage is 
reached, the molten steel is covered with a hot slag yielding conductive vapors that help 
to stabilize the arc current and the load on the supply system. 

Empirical data collected by the utilities show a great variance in single-phase load swings 
among furnaces [7]. Abrupt changes in the circuit impedance cause these load swings. 
The impedance of the arc part of the circuit will vary for a number of reasons such as 
those listed below, 

1. The metal in the furnace may melt away from the electrode or it may be forced away 
from the electrode by the inductive effect of the arc current, thus abruptly increasing 
the impedance of the arc. 

2. The electrode oxidizes away and ionizes. This increases the arc impedance slowly. 

3. The physical and chemical make-up of one charge is different from the next. This is 
known to have a great effect upon the behavior of the arc. 

4. The arrangement of the charge materials in the fiimace has an effect upon the 
behavior of the arc. 

5. The relative amount of air blown upon the arc and its moisture content will have a 
greater or lesser chilling effect to increase or decrease the impedance of the arc. 

Thus it can easily be seen that electric arc furnace is a large and potentially severe load 
on the electrical supply network. Owing to the nature of the electric arc phenomena, and 
depending on the supply capacity, the size of the furnace and the mode of furnace 
operation, the arc furnace can result in serious levels of electrical disturbance on the 
network. Estimation and control of the disturbances is very much a system problem, and 
at the design stage of a new installation, input from the electric utility, the furnace 
manufacturer and the plant engineering staff will be required, if the level of disturbances 
is to be kept within limits. 
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1.3 Methods of arc fiirnace compensation 

Conventionally capacitors were used in the electrical system that acts analogous to an 
accumulator in a hydraulic system. They will draw and store the excess energy from the 
system whenever the force producing the energy exceeds normal and return this energy to 
the system whenever the force fells below normal. Thus capacitors are a means used to 
help stabilize the voltage and current in line. 

A synchronous condenser may be designed to run on active or reactive energy from the 
line and produce active or reactive energy in the form of mechanical energy or electrical 
energy. The mechanical energy may be drawn off the shaft of the rotating machine by the 
attachment of a gear. The electrical energy, as active watt-hours or reactive var-hours is 
fed back into the line. 

The arc furnace is a highly reactive load because of the inherent and very often added 
reactance in the circuit to stabilize the arc. The arc fiimace draws energy from the line, 
converts only a part of it into useful heat, and the rest is fed back into the line as reactive 
energy. It does no good to the line, and a synchronous condenser is used to convert it 
back to useful energy again. By converting the reactive energy to active energy again, the 
synchronous condenser converts a low power fector load to a high power factor on the 
line and its flywheel effect absorbs fluctuations in the line voltage and current caused by 
the arc furnace. 

However in recent years, static VAR compensators using thyristors have been built which 
are fast superseding the conventional methods for arc furnace compensation because of 
their being superior in respect of response, maintenance and efficiency [7]. They can 
compensate for the fundamental reactive as well as, to a limited extent, the harmonic 
content in the supply current. 
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1.4 Problem formulation 


It is seen that an arc fiimace behaves quite randomly in the sense that the arc impedance 
changes abruptly. The real and reactive powers drawn by the furnace also fluctuate 
randomly and often severely thus putting considerable stress on the system. This leads to 
a poor voltage profile at the bus supplying the arc furnace installation and also creates 
flickering of incandescent lamps and other disturbances to the other loads being supplied 
jOrom the same bus. 

To overcome these problems, different compensating schemes have been employed, the 
one that is being used at most of the installations is compensating with the help of 
S.V.C’s . Although they are far better than the synchronous condensers that were 
conventionally used but still they are not able to compensate for the problem of flicker to 
reasonable extent, as will be shown in the following chapters . 

The objectives of this study are, 

> Formulate an appropriate model for an arc so as to simulate a furnace installation. 

> fivaluate the performance of the arc furnace, 

1 . Without any compensation. 

2. With presently used S.V.C scheme. 

> Investigate the affects of compensating the arc furnace instantaneously with the help 
of a controlled current source on the reduction of flicker and other problems. 
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1.5 Thesis outline 

This study has been performed keeping the fore mentioned objectives in mind. The entire 
study revolves around how to simulate an arc furnace and then compensate for the flicker 
and other problems. 

The thesis is organized as follows, 

Chapter 2 deals with different models of the arc and the one that has been structured so as 
to cover different aspects of arc furnace operation used for the purpose of this study. This 
chapter also gives the benefits of the proposed model as compared to the conventionally 
used models of an arc. 

Chapter 3 gives an idea of the kind of load an arc furnace is, looking from the source 
point of view. An arc furnace installation has been simulated and the behavior of the arc 
has been studied. The chapter gives an idea of the kind of harmonics generated and the 
reactive power required by the arc furnace. Flicker has also been defined and evaluated 
for the bus supplying the furnace installation. 

Chapter 4 discusses the performance of the furnace installation with a S.V.C connected in 
parallel for the compensation of reactive power. It is shown that although we have a 
reduction in the average reactive power required by the furnace but still the harmonics in 
the source current are not suppressed to the required extent. What is most important is the 
level of flicker after compensation. It is shown that we do not have appreciable reduction 
in the level of flicker at the critical bus (the one supplying the fiirnace). 

Chapter 5 revolves around the theory of instantaneous compensation proposed by Lai 
[10] and the modification carried out to take care of the harmonic powers as well. The 
method used to evaluate the reference current for compensation takes into account the 
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fundamental as well harmonic powers drawn by the load thus giving far better results for 
rapidly varying loads such as arc furnaces. 

Chapter 6 investigates the effects of compensating an arc furnace with the help of, first an 
ideal current source and then with a current source realized by a 21evel, 3phase voltage 
source inverter. It is shown that by using the scheme developed in chapter 5 and an ideal 
current source we can in fact improve the performance of the furnace excellently. We’ll 
see that the harmonics in the source current are reduced drastically and the level of flicker 
comes down by as much as 80-85%. Although the performance deteriorates slightly when 
compensating with a 21evel inverter but still the results are far better than those obtained 
with S.V.C. 

On an overall the study brings out the fantastic improvement in the performance of arc 
furnaces, if compensated instantaneously with voltage source inverters as compared to 
Ihe presently used compensallon Hchemes using S.V.C's.. As concluded In chupter 7, we 
can see that still there is scope for improvement upon the results obtained in this study 
and some of the possible alternatives have been listed as the scope for fiiture work. 
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Chapter 2 


MODELLING OF AN ARC 


2.1 Introduction 


One of the major difficulties in simulating an arc furnace installation is the availability of 
a proper and concise model of an arc. Although arc has been studied to a great extent but 
still a satisfactory model is lacking basically because of a highly non-linear relationship 
between the arc voltage and ciuxent and also due to the unpredictable behavior of the arc 
current during the initial period of the melt. 

Thougl^ we have in literature [9] quite a few models of arc but their main drawbacks are 
that either they are valid for relatively small current area (Circuit breakers, Welding 
equipment) or they involve thermodynamics to a great extent thus requiring elaborate 
knowledge of thermodynamics and material science. 
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We present here three models of an arc that does not require much background of thermal 
physics. 

1 . Electro Magnetic Transient Program model. 

2. Arc current simulation model. 

3 . V ariable resistance model. 


2.2 Electro Magnetic Transient Program Model 


A complete three phase model can be implemented using E.M.T.P [14]. The model of an 
arc is given by a modified version of the Mayr equation, the Cassie-Mayr formula 



Rare 

~~er 


V^arc 


V^o 


where, 

> Rare is the arc resistance. 

> ^ is the arc time constant. 

> Vo is the arc static voltage depending upon the arc length. 

> Varc is the arc voltage. 

Such an implementation can be done using E.M.T.P. Transient analysis of control 
systems module and the use of controlled voltage or current sources. 


2.3 Arc Current Simulation Model 


One way to simulate an arc furnace is to model its current waveshape [3]. A number of 
studies and measurements have been carried out into the arc furnace parameters and their 
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harmonic spectrum based on which a computer model for an arc current can be 
constructed. A typical arc current model with predefined magnitudes of different 
harmonics and a random variation about the same has been given below. 

^arc (^) = sin(wt) p.u 

^arc(0 = Ia(t) + Ir(t) P-^ 

la ~ ( 0.5 X I + 0.75 ) sin ( w t ) p.u 

(t) = -(X 2 +0.5)cos(wt) + (X 2 -0.5) 

+ 0.22 {(X 3 -0.5)sin(2wt) + (X 4 -0.5)cos(2wt)} 

+ 0.26 {(X 5 -0.5)sin(3 wt) + (X^ - 0.5 ) cos ( 3 wt) } 

+ 0.10 {(X 7 -0.5)sin(4wt) + (Xg -0.5)cos(4wt)} 

+ 0.22 { ( X 9 - 0.5 ) sin ( 5 w t ) + ( Xjq - 0.5 ) cos ( 5 w t ) } 

+ 0.04 { (Xj! -0.5)sin(6wt) + (Xj 2 - 0.5 ) cos ( 6 w t ) } 

+ 0.12 {(Xi 3 -0.5)sm(7wt) + (Xi 4 -0.5)cos(7wt) } 

+ 0.04{(Xi 5 -0.5)sin(9wt) + (X 16 - 0.5)cos(9wt)} p.u. 


where, 

is the arc voltage. 

I is the arc current. 

Xi ” Xjg are random numbers distributed uniformly between 0 and 1. 
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2.4 Variabie Resistance Model 

It has been shown in literature [11] that an arc essentially consists of a fixed reactance in 
series with a variable resistance. Hence we can also model an arc if we can vary a 
resistance randomly about a mean operating point in such a way that an arc current can be 
simulated. 

Arc current is generally seen to have large fluctuations particularly during the initial 
period of the melt. Hence we simulate the arc by switching in and out different fixed 
resistances about a mean value and further randomly select the time period of their being 
in the circuit in such a way so as to best simulate the actual arc current and its harmonics. 


Lf 



Figure 2.1 : Arc furnace model 

Above figure shows how a combination of resistors and switches can be connected to 
simulate an arc furnace. Only one switch is on at a time and that switch is selected 
randomly 

Table 2.1 shows the values of resistors used for simulating the arc, the switch positions 
and the no. of time steps for which that state will remain before another resistor is 
selected . 
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SI 

S2 

S3 

S4 

S5 

Rare 

Time 

R1 =3 

ON 

OFF 

OFF 

OFF 

OFF 

3 

X*100 

R2 = 8 

OFF 

ON 

OFF 

OFF 


8 

X*300 

R3 = 10 

OFF 

OFF 

ON 

OFF 

jS^jBjl 



R4=15 

OFF 

OFF 

OFF 

ON 


15 


R5 = 20 

OFF 

OFF 

OFF 

OFF 

ON 

20 



Table 2.1 


In the above table ‘ X ‘ , denotes a random number uniformly distributed between 0 and 

1 . 

This method has following advantages over others and hence is being used for the 

purpose of simulating arc. 

> There is enough unpredictability in the value of the arc resistance to come in the 
following time as well as the duration for which it will remain thus simulating to a 
great extent the dynamic behavior of the arc . 

> The arc current spectrum can be designed so as to match that of an actual arc furnace 
by appropriate weightages given to different resistance values. 

> The power drawn by the arc can be varied by having a different set of resistance 
values . 

> We have time intervals during which the arc has a negative resistance characteristic 
due to sudden variation in the resistance of the arc thus simulating the actual behavior 
of an arc. 
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Chapter 3 


ARC FURNACE AS A LOAD ON THE UTILITY 


3.1 Introduction 


An electric arc furnace, as such, is one of the most troublesome loads to be supplied by 
the utilities, but looking from commercial point of view they are large revenue generators 
and hence form important loads for the utility. The basic problem with the arc furnace is 
the randomness associated with the arc impedance owing to which the arc current tends 
to fluctuate considerably from cycle to cycle. The only check to these rapid fluctuations 
is the reactance put in series with the arc, thus limiting the rate of change of current and 
also limiting the short circuit currents which is a common incident during the early 
phases of the melt [11]. The reactance forms an essential part of the system but on the 
other hand consumes large amounts of reactive power thus putting considerable stress on 
the system and limiting the real power flow to the arc furnace. 



14 


Arc furnace installations generally have large transformers that have high leakage 
inductance. As the arc draws large and fluctuating currents the voltage profile at the bus 
feeding the arc furnace deteriorates and the voltage fluctuates giving rise to the 
phenomena of flicker which is defined in the next section. 

This chapter highlights various problems and the phenomena of flicker associated with 
the arc furnace if operated vvithout any preventive and curative methods. The variations 
of different parameters such as, current, real and reactive power, voltage fluctuations and- 
their harmonic spectrum are also studied. 


3.2 Flicker 


Flicker is defined as the low frequency (0.1 Hz to 20 Hz) modulation of the supply 
network voltage. The modulation results in characteristic flickering of incandescent 
lamps which can be most annoying when the flicker is in the 1 to 10 Hz range. Very 
small variations are enough to induce lightning disturbances unbearable to the human 
eye. For a standard 60W lamp, the disturbance becomes perceptible for a voltage 
variation frequency of lOHz and relative magnitude of 0.3%. 

Whenever an arc furnace is connected to a supply having a reasonable short circuit level, 
a proper assessment of the resulting flicker level is probably the most important aspect of 
the disturbance question to be addressed. 

A common indicator of the level of flicker is obtained on a scale called A V jq [3], which 
is generally used in Japan. A Vjq is the rms value of the voltage fluctuations weighed by 
the frequency characteristic of the human eye’s sensitivity for flicker and is calculated by 
the following equation. 
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1/2 


^^10 = 


3.1 


where , A V „ : The frequency component at frequency f „ of the voltage . 

a„ : Flicker sensitivity coefficient at frequency f„ of the human eyes as 
given in [4]. 


3.3 Arc furnace system outline 

The single line diagram of the system consisting of the step down transformer the arc 
furnace transformer and the arc furnace is as shown in figure 3.1. 


R. L. 



Vs 


Tr. 1 



Tr. 2 


(■ 


r' 

Arc r " 
Furnace 


Figure 3.1 : Single phase equivalent of the arc furnace installation 

Here , 

Vs= 154kV,3phase,60Hz 

Rs = 0.254 Q , Ls = 58.1 //H (System Impedance). 

Trl=154kV / 33kV , 91MVA , X = 13 % . ( Step Down Transformer ) . 

Tr2= 33kV / 812V , 75MVA , X = 7.5 % ( Arc Furnace Transformer ) . 

Laf= 30 mH ( Current Limiting Reactor ) . 
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The 33kV bus between the step down transformer and the furnace transformer is being 
considered as the critical bus whose voltage fluctuations causes flicker at other 
installations on the same bus 


3.4 Single phase model and system equations 

Due to the randomness associated with the arc the arc furnace essentially behaves as three 
single phase loads particularly during the early phases of the melt, hence it can be said 
that all the phases are relatively independent of each other. The single phase equivalent of 
the system is shown in figure 3,2. 


Is 



Figure 3.2 : Single phase equivalent of the arc furnace system 

The values of all the parameters can be obtained from the data on the last page. I’he value 
of is chosen such that the short circuit current is limited to around 3 kA. The value 

of R is variable and varies between 3 to 20 ohms in a fashion described in the previous 
chapter. 
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The differential equations governing the system are given by, 


~ i^s + + Al + ^2 + ^qf ) 


dt 


3.2 


and 




dis 

~dt 


3.3 


The differential equations are now converted to linear state equations of the type , 


and 


X=AX + BU 

Y=CX + DU 


3.4 


3.5 


where , 

X = Is , Y = Vp , U = Vs. 


{Rs + 

1^5 +'^t2 +^af) 


B = 


+ ^/2 ) 


r 


C^-Rt 


1 + 


+^tl +^t2 +^a/) 


( 


D 


(^s +'^rl +'^/2 +^qf) 


The system can now be solved on Matlab by converting it into a discrete model in the 
time domain . 



18 


3.5 Variations of different parameters during initial period of the melt . 

The system is solved for a period of 1.5 sec. so as to obtain harmonic spectrum upto a 
frequency of IHz required to evaluate flicker on the critical bus , the different waveforms 
and their analysis is given below . 


3.5.1 The Arc current and its harmonic spectrum . 

The Arc current fluctuates considerably and has a large harmonic content as shown in the 
figure 3.3 .It is quite clear from the figure that the current has some low frequency 
random modulations which are the main cause of flicker .The rms value of the current 
also changes randomly leading to a deteriorating voltage profile at the critical bus . 



Figure 3.3 : Arc current waveform 








Current ( *0.1 A ) , Voltage ( V ) 
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X 10^* 



Figure 3.4 : Load current and source voltage 

We can see from the figure 3.4 that the source current lags the voltage by a large amount 
and hence there is strong need of compensation of some type to reduce that level.The 
spectrum of source current shown in figure 3.5 also highlights the strong need of 
compensation . 



Figure 3.5 ; Load current spectrum 
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Figure 3.6 indicates that the real and reactive power drawn by the load fluctuates 
tremendously thus highlighting the importance of the requirement of some form of 
compensation or else the supply system may suffer from disturbances . 

3.5.3 Voltage at the critical bus . 

The following plots shows the variation of the voltage at the critical bus with respect to 
the supply system voltage . 


X 10^ 



Figure 3.7 : Source and critical bus voltage 


It can be seen from the figure that the phase voltage can drop by as much as 5 % of the 
rated supply voltage .Further it is clear from the next plot that the voltage fluctuations 
have a high harmonic content at low frequencies particularly below 30 Hz , giving rise to 
flicker . 
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3.5.4 Evaluation of flicker on the critical bu» ■ 

Flicker has been already defined in the previous section. Taking that definition and the 
above spectrum for voltage fluctuations we have 

fao 

where, 

= Sensitivity coefficient of human eye at frequency fh . 

AF„ = Harmonic amplitude at frequency fri . 

The flicker level level for this uncompensated case turns out to be 58.5 . 


3.6 Conclusions 


It is seen that the arc furnace indeed behaves in a very random fashion. The arc current 
contains a large amount of harmonics thus leading to a poor voltage profile at the critical 
bus. The power factor is also very poor and the level of flicker is very high. 

Thus we can conclude that an arc furnace just cannot be operated without any 
compensating elements. It has also been shown that the furnace parameters vary 
randomly and at times very quickly, hence we need a compensation scheme that is fast 
enough to compensate the arc fiirnace effectively and reduce the flicker to acceptable 
levels. 
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Chapter 4 


COMPENSATION OF ARC FURNACE WITH S.V.C 


4.1 Introduction 


The capacity of arc furnaces have risen at a sharp rate in the recent years, as a result, the 
problem of flicker has become quite acute. Traditionally, rotating synchronous 
condensers were used in order to compensate for the reactive power and reduce flicker. 

In recent years however, static V. A.R. compensators using thyristors have superceded the 
synchronous condenser for the following reasons, 

1 Fast response (within a lialf cycle). 

2 Low maintenance costs. 

3 Superior efficiency and lower losses. 
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V.A.R compensation using thyristors have the following benefits besides flicker 
reduction. 

1 Improves the overall power factor. 

2 Reduces the harmonic currents. 

3 Reduces losses by eliminating the reactive current. 

This chapter investigates the effects of a shunt installed S.V.C., comprising of a fixed 
capacitor and thyristor switched reactor, on the operation and characteristics of arc 
furnace installatioBu 


4.2 Basic performance of a S.V.C. 


The basic idea of S.V.C. is to control the reactive component of the source current, 
through compensating the load, by controlling the current through a reactor with the help 
of thyristor switches. Two anti-parallel thyristors Thl and Th2 are connected in series 
with an inductor L, as shown in the figure 4.1. We also have a capacitor in parallel with 
this reactor branch so as to make it work as an effective capacitor. 


icx>m ^ H 



Figure 4.1 ; Principle of operation of a S V.C 



Current ( A ) , Voltage ( V ) 
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Thyristors Thl and Th2 are made to conduct by means of a trigger signal, the duration 
and peak of current is regulated by the variable angle of delay of the trigger signal, a , 
with respect to the peak of the phase voltage. This is shown in the figure 4.2, 


X 10^^ 



Figure 4.2 : Thyristor current for different firing angles and source voltage 

Hence we can say that with the help of the thyristor switch we can control the reactor 
current and in effect control the reactive power drawn by the inductor. 

The reactor current at a firing angle delay of a with respect to the voltage peak is given 
by the following equation, 

ij(t)= if a^W't^{n;-a) 4.1 

I 

where , ; 


Iiq = — — = peak value of reactor current at zero firing angle . 
w 'L 

Vq = peak value of the supply voltage . 




27 


The fimdamental component of the reactor current at firing angle a is given by the 
following expression. 




4.2 


The firing angle at rated reactor current is called . The selection of depends upon 
the tradeoff between the flicker compensating effect and harmonics generated by the 

reactor itself. For the purpose of this study has been taken as 30^ this implies that 

the firing angle can now vary from this value to 90^ . 

Using a fixed capacitor in parallel with the thyristor switched reactor makes the complete 
unit work as an effective variable reactive power source. The total compensator current is 
the resultant of the capacitor current and the reactor current and is given by, 

4.3 

where , 


^cW=-^co sin(>v-0 

where I • tV • C 

and is given by the previous equation . 

The fundamental component of the total compensator current is given by. 


^comf 


(a)=Ic 




+ 





4.4 
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4.3 Measurement and control 


In this type of V.A.R compensators, it is important to measure the reactive power 
variations of the arc furnaces precisely and control the thyristor qurrent as rapidly as 
possible to accomplish good compensation. 

The principle of measuring and control state that the active current ip of the arc furnaces 

be measured by multiplying the instantaneous value of i at the peak of the voltage by a 
unit sinusoidal wave in phase with the voltage and the reactive component can be 
detected continuously by subtracting ip from / i.e. iq = (/ — ip) • 

Another way to measure the required fundamental reactive current Iq is that we continue 

evaluating the average reactive power demanded by the load for the last half cycle by the 
moving average technique and estimate the required reactive current by dividing the 
reactive power required by the load with the rms voltage at the bus . 


When the furnace current changes abruptly the measuring signal TQ’, equal to the 
fundamental component of the reactive part of furnace current for the last half cycle, also 
follows the changes immediately. The firing instants for the next half cycle are generated 
by evaluating the crosspomt of TQ’ and the function ‘F’ wWch represents the relation of 
the fundamental compensating current and the firing angle as shown below in figure 4.3, 
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Figure 4.3 ; Comparing IQ and Icompf obtain the firing angle 

Thus it is clear that we can compensate for the fundamental reactive current of the arc 
furnace for each half cycle in the next half cycle. 


4.4 System outline and state equations 


The system that has been considered for compensation in none other than the one 
considered in chapter 3 except for the addition of an S.V.C on the 33kV bus. The figure 


Rs L, 




Figure 4.4 : Single line diagram of the arc furnace installation 
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4.4 shows the single line diagram for the arc furnace installation being considered for 
compensation. 

The single phase equivalent of the system along with the values of all the parameters is 
given in the figure 4.5 shown below. 



Figure 4.5 : Single phase equivalent of the arc fiimace system 


The values of S.V.C parameters are as follows, 

Cf = 0.1 mF Lf = 5 mH L = 50mH 

The differential equations governing the system are as follows, 

Case 1 : If any of the thyristors is conducting , 

K = i^s +^n)’ h + i^s + ^t\ )-^ + (^) * V + Lqf )-^ 4.5 

)• h + (^5 + ^tl )~^ + L 


4.6 
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^ s = )• L + i^s + ) 


di 


dt 


' +yc+i^f 


dt 


4.7 



dt 


4.8 


Case 2 : When none of the thyristors is conducting , 

All of the above equations are valid except the second. But then the number of state 
variables also reduces to three because of elimination of Ix. Hence in that case also the 
system can be solved with three equations at hand. 


The above equations can be solved by converting them to state equations of the type , 

X = AX + BU 

4.9 

Y = CX + DU 

where X is the state vector with four state variables and U is the input function that is the 
system voltage , A , B , C , D are system matrices which are given below , 


X = 


ri/ \ 


V 






32 


further , 

if none of the thyristors is conducting , 


A = 


0 


0 


1 


]\>^Y ^^2. ^^3 

(jLj'My — l) (Z/2’A/^2 "h jR "h jRj’ ) Lj'AA'^ 


0 


0 


Lf 

0 


0 

0 

0 

0 



0 


M5 . 

B = 





0 


and if any of the thyristor is conducting , 


0 

0 0 0 ' 


0 

Ml 

M2 M2 M4 


M5 

{LjMi -l) 

^Lj' At 2 R Rj''^ Lj'At 2 LyM[ ^ 

B = 

LtM^ 

h 

Lf Lf Lf 


w 

LjpAd j 

(Z/j’A/2 "b Jf? "b Z?j’ ) Z/j’A/j Lj'M/^ 



L 

L L L . 


_ L _ 


and for both of the above cases , 


Lj'M\ 

Lj'Kt^, "h jR "h Rj' 
X j'A /^3 
Z!/ j’ A jT 4 


D = [LrM^] 
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where , 



Thyristors not conducting 

Thyristors conducting 

D 

r r f 1 1 0 

/ / - 1 - 4 

f 

L^Lt 

\ 

1 1 1 0 

JU V JL/'T 1 T 

\h irj 

1 I 1 

Lf Lg Lj Lj 


Ls 


Rs 


LfD) 

[h-D] 

M 2 

1 „ 

1 "i” 

D ^ 

V. 

^ L ^ 

— + 1 (i? + /27<)l 

J J 



M3 

D 

Rs 

D 

M4 

0 

Rs 

D 

Ms 

1 

D 

1 

D 


Here in the above expressions, 

^5 = ^5 + ^/1 ^5 = + ^ t \ 

Rt = Rt2 


- ^/2 + Lqf 
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Update supply voltages 

sm{wt) 

Vj, = V„sm{wt-m) 

Fc=FoSin(w/-f-120) 


Update supply real and reactive power { & q^) 

Update load real and reactive power (. Pi &qi) 


J = 0; 

Update arc resistance , , and 

system matrices 
■^o ’ 

Ai j B2 i C2 ’ ^2 
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4.5.4 Flowchart for the function UpdateAloha 







Arc current , Source current 
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4.6 Performance of arc fiimacc with S.V.C in parallel 

The rest of this chapter deals with evaluating the performance of the arc furnace with a 
S.V.C in parallel. 

4.6.1 Arc current and source current waveform and sncctnini 

The waveform shown below indicates that we do have significant reduction in the mean 
r.m.s value of the current being supplied by the source basically because of compensation 
of the average reactive current required by the arc furnace. It can also be seen that we do 
not have very good compensation as the source current, though less in magnitude, still 
varies quite a lot and has a large harmonic spectrum as indicated in the figure 4.6 given 
below. 
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0 



Figure 4.7 : Load and source current along with source voltage 
It can be clearly seen from the figure 4.7 given above that the supply current, though not 
sinusoidal indeed is almost in phase or leads slightly owing to the capacitor connected in 
parallel. This is probably the best we can do with a S.V.C because it cannot be controlled 
‘ in less than a half cycle and hence is inherently incapable of compensating ideally such a 
rapidly changing load. 
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Figure 4.8b : Source current spectrum 

From the figures 4.8a and 4.8b given above we can see that we do have some limited 
improvement in the spectrum at frequencies other than 60 Hz. We are particularly 
worried at lower frequencies because they are the main cause of flicker problem. Hence it 
can be asserted that a S.V.C is a bit slow to efficiently compensate an are furnace. 


4.6.2 Real and Reactive power 

The figures 4.9 and 4.10 given below indicate the signific^t reduction in the demand of 
reactive power from the supply as a result of compensation with a S.V.C .It can also be 
seen that the fast variations in the reactive power demand are not suppressed to a great 
extent although on an average basis the demand has reduced drastically . 

Thus it can be said that though a S.V.C might be a good compensator for other loads but 
it is definitely not so with arc furnaces. The variations for arc furnaces are wild and 
random, and hence they are far too fast for the S.V.C to respond. 





Load and source real p 



0.4 0.6 0.8 1 1.2 1.4 

Time ( Sec. ) 


Figure 4.10 ; Load and source real power 
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As for the real power, it also varies randomly along with the arc resistance. As we have 
compensated only for the fondamental reactive power hence all the real and reactive 
harmonic power is still derived from the source hence resulting in such a waveform. 


4.6.3 Voltage and its spectrum at the critical bus 

As shown in the previous chapter, the voltage at the critical bus fluctuates whenever the 
arc fiirnace is in operation, thus resulting in flicker, ^^en an S.V.C is connected in 
parallel to the furnace the voltage profile improves significantly though the harmonics 
present in the voltage spectrum are not reduced to a great extent. 



0.62 0.64 0.66 0.68 0.7 0.72 0.74 0.76 0.78 


Time ( Sec. ) 


Figure 4.11 ; Source and critical bus voltage waveforms 


Hamionic amplitude Voltafle at critical bus ( V ) 
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The figure 4.12 given below shows the variation in the r.m.s value of the voltage at the 
critical bus and the corresponding voltage spectrum at the critical bus with a S.V.C 
connected in shunt with the furnace is given in figure 4.13. 


X 10^ 




Figure 4.13 : Harmonic spectrum of the critical bus voltage 
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By comparing this figure with the corresponding figure in the previous chapter, it can be 
seen that we do have some improvement in the spectrum but not very significant and 
hence we can say that the flicker will not be reduced to a great extent as can be seen in 
the next section . 


4.6.4 Flicker at the critical bus 


Flicker is calculated at this bus in the same way as described in the previous chapter. As 
has already been seen from the previous figure that the voltage spectrum does contain 
quite a large amount of harmonics we do not expect a great improvement in the level of 
flicker at this bus. 

Using the same definition of flicker level i.e. AF^q we achieve a flicker level 

AFio =50.12 

Comparing this with the level obtained in the previous chapter we do see an improvement 
of around 15%, but it is limited primarily because of the slow speed of the S.V.C as 
compared to the fluctuations in the arc furnace parameters. 


4.7 Conclusions 


Having seen the operation of the arc furnace first without any compensation and then 
with a S.V.C connected in parallel we can conclude the following, 

> The average reactive power required by the arc fiimace can be compensated to a great 
extent by a S.V.C so that the load seems to be a resistive varying load if talked on an 
average basis. 


46 


> The source current though reduces in magnitude owing to the compensation of the 
fundamental part of the reactive current but still it has a lot of harmonics because of 
the discontinuous operation of the thyristor. 

> The real power consumed by the load also shows large harmonics indicating the 
presence of large amount of real harmonic power which cannot be compensated with 
the help of a S.V.C and hence we ought to have a different method of compensating 
the reactive current which is fast and whose response is, ideally speaking 
instantaneous . 
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Chapter 5 


THEORY OF INSTANTANEOUS COMPENSATION 


5.1 Introduction 

For single phase power systems with sinusoidal voltages and currents, quantities such as 
active power, reactive power, active current, reactive current, power factor etc. are based 
on the average concept. 

Many authors have attempted to redefine these quantities to deal with 3 phase systems 
with unbalanced and distorted currents. Among them Akagi et. a/ [ 1 ] have introduced an 

y 

j 

interesting concept of instantaneous active and reactive power. This method provides an 
effective means to compensate for the instantaneous components of reactive power for 3 
phase systems, but the theory does not work properly in the presence of zero sequence 
components as pointed out in [10]. As the arc furnace does contain zero sequence 
currents, we have for the purpose of this study used the theory proposed by Lai [10], 
which is explained in the next section. 
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5.2 Instantaneous real and reactive powers 

For three phase systems, both the currents and voltages can be treated as vectors with 3 
components one for each phase v.i.z. 



fi ^ 

*a 


fv "i 
(1 





II 

'Jc ^ 

; V = 

n 


5.1 


where , , ijj , are the instantaneous values of phase currents . 

Va,Vb> instantaneous values of phase voltages . 

Now the theory proposed by Lai defines the instantaneous real and reactive powers as 
follows 


p = V 'i 



5.2 


Applying simple vector algebra rules we have 


where. 


P = ■‘'a • *0 + ^6 ' h + ' h 


5.3 


q 


{ql+ql+qc 



5.4 


h h 



qc 


V. 


h 


qa = 



49 


Taking these as the fundamental definitions it can easily be shown that the instantaneous 
current vector can be decomposed into two parts one being the instantaneous active 
component and the other being the instantaneous reactive component as follows 




where , 




r p \ 


V 


it can also be shown that 


q XV 
(v-v) 


5.5 


5.6 


where, 


s = v>i ; X = — 

s 


(i " 

^ap 

% 

f/ ^ 

Hp 

: ?,= 

hq 



\Jcq) 

H|v1l 

and i = 

_*.! 

1 


5.7 


s = instantaneous apparent power . 
A= instantaneous power factor . 


Thus it can be argued that the load can be compensated if we can somehow compensate 
the reactive part of the current drawn by the load. 
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It turns out that the proposed theory works well in the case of balanced voltages but not 
in case of unbalanced voltages which is precisely the case with the arc furnace. To 
overcome this problem we consider the scheme proposed by Watanabae [13]. 


The active and reactive power consumed by the load in general varies about some mean 
value and hence both of them can be decomposed into their respective mean and 
oscillating components i.e. 


p=p+p ; q=q+q 


5.8 


where, 

p, q are the average values of the instantaneous real and reactive power respectively. 
p, q are the oscillating components of the same . 

The concept of instantaneous compensation involves that the source should supply only 
the average real power and everything else should be compensated for as depicted in the 
figure 5.1 given below. 


M — ► 



Figure 5.1 : Process of compensation of different power components 

As is obvious from the definitions, these power components cannot be achieved directly. 
The process of evaluating them is being given in the next section. 
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5.3 Moving average technique 


In order to estimate the average and alternating components of the real and reactive 
powers we either have to use a low pass filter or consider a moving window that 
considers only the latest values of the parameters. The second method has been employed 
for the study as it is faster and moreover does not effect the compensation characteristics 
because of instantaneous response unlike other low pass filters. 

The averaging of p and q is performed over a period of one half cycle of the supply 
voltage as both p and q have fundamental components at twice the supply frequency. The 
moving window considered has the latest sampled data of instantaneous real and reactive 
powers for the last half cycle. These data are used to find the average of p and q i.e 
p • Once we have the values of p Sc q we can subtract these from the 

instantaneous values of p Sc q respectively to obtain p Sc q .The process has been 
explained below. 


Suppose we have N samples in the window and suppose the parameters for the 
sampling instant are being evaluated then, 


= < 






[V-^ )k={n-N+\) 


whenever n'^N. 5.9 


p = p-p , 5.10 

As the time increments so does our moving window and so does the power components. 
Thus at any time we have all the power components derived from the data for the past 
half cycle and that is the best we could do. 

The next section deals with the process of generating the required reference currents so 
that we may compensate the entire reactive and the harmonic real power. 

:EmiiKL tmup 
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5.4 Reference current generating 

Suppose that the voltage vector at the point of compensation is given by V; and the load 
current to be compensated is given by Ij then as pointed out earlier. 


Pi =Vi ■ ii and qi=ViX Tj 

further, 

(v/-V/) (v/-V/) 


On simplification we have. 


V 



5.11 


5.12 


According to Lai if we can compensate this current we will have a virtually balanced 
3phase load thus. 


5.13 

(v, -v,) 

1 

but as the arc furnace is typically one of the worst kind of loads the system is being 
subjected to and as has been shown earlier , the reactive as well as the real power being 
consumed varies wildly , the £*ove mentioned theory , does not give very good result 
basically because the method relies on pompensating the reactive current term and hence 
the arc furnace continues to taking the rapidly oscillating real power from the supply thus 
containing harmonics in the supply current . 
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Now we compensate for the harmonic real power by appropriately modifying the value of 
reference current required as follows 


qiXV 


hef 


Vi ‘V; 




5.14 


where Pi is the alternating component of the load real power . 

Thus now the source has to supply only the average real power and hence we have a 
considerably better compensating scheme although it has some problems due to the 
) random nature of the arc current but still this scheme turns out to be far belter than the 
methods used presently for the purpose of arc furnace compensation . 

On simplification the above expression becomes 


^ 7 V 


5.15 


where, 


Pi is the average value of the load real power and is equal to pj -pi and is calculated 
by the moving average method as explained earlier. 


Now if we can somehow truck this rcfcronco current, the entire reuelive and llic hunnoiuc 
real power will be compensated. Although this does not guarantee perfect compensation 
because the reference current has been generated with the help of data for the last half 


cycle and the furnace parameters may vary at the time of compensation. 
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S.5 Conclusions 


Thus it can be concluded that the load currents and powers can be decomposed into their 
respective average and oscillating components namely, 

PhPi & q^qi 

It can be said that Pi &qi contribute towards the current harmonics present in the 
supply current . In normal cases the contribution of is far greater than that of pi and 
hence theory proposed by Lai works very well but in the case of arc furnace pi is far 
greater than that for other equivalent loads and hence it is absolutely essential to 
compensate for g'/ as well as pi . 

The implementation of this theory now involves that these reference currents be 
continuously tracked and compensated for with the help of suitable current source an4 
control scheme such as bang-bang hysteresis control. 
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Chapter 6 


COMPENSATION WITH A CURRENT SOURCE 


6.1 Introduction 


Several different S.V.C configurations are presently used for compensation of arc 
furnaces. All the conventional S.V.C’s employ thyristors to vary the effective value of 
passive reactive elements connected in parallel with the load. The design of these S.V.C’s 
is such that the compensation of the load occurs discontinuously in steps of one half cycle 
or more of the power line frequency. Conceptually S.V.C’s are intended to compensate 
only the fundamental reactive current of the arc furnace and hence the latter continues to 
take the harmonic real and reactive powers from the source, which is the main cause of 
flicker. 

Tuned filters are also used in conjunction with the S.V.C’s if harmonics in the load 
current are high .The inherent delay of these S.V.C’s may be too great for compensation 
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of loads such as arc furnace where the behavior of the arc is a random phenomena. It has 
already been justified in the previous chapters that the conventional methods of 
compensating the arc furnace with the help of S.V.C’s falls short of the requirements 
necessary for flicker management. This is primarily because of the fact that the arc 
" current changes rapidly and so does the required compensating current. Thus it reasonable 
to investigate the compensation effects based on the theory proposed in the previous 
chapter, by using inverters as current sources to compensate the arc furnaces 
instantaneously. 

This chapter deals with compensation of arc furnace system instantaneously with the help 
of current sources. The current required for con^nsation is derived instantaneously and 
the current source tracks this reference current in order to accomplish compensation. This 
method has some distinct advantages over the conventional methods using the S.V.C’s 
such as, 

> The response is very fast compared to conventional methods. 

> Not only the fundamental but also the harmonic reactive currents can be 
compensated. 

> By properly generating the reference currents we can also compensate the harmonic 
real power so that the arc furnace virtually behaves as a balanced resistive load. 

The first part of this chapter deals with compensation of arc furnace with an ideal current 
source and the second part deals with compensation with the help of a current source 
realized with a 3phase 21evel inverter. 


6.2 System conlRguration and state equations 

The system considered for compensation is again the same with the only difference that 
now we have a current source in parallel with the furnace as shown in the figure 6.1, 


57 



In the above figure is the current supplied by the compensator and follows a 
generated reference current as explained in the next section. 

The equations governing the arc furnace system shown above are, 

dl 

Vp=(Ra*R)-I^*{L,2^L^)~ ' «-2 


where, 

is the co mp ensator current that follows the reference current which in turn is 
generated as explained in the previous chapter . 


These set of differential equations are again converted to state equations of the type. 



X = AX + BU 

and 

Y = CX + DU 
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6.4 


where, 



In the above expression is a function that represents, 

> The max. rate of change of current through the current source in case the compensator 
is assumed ideal . 

> If the compensator is assumed to be of 3phase 2 level type, this represents the 
instantaneous voltage difference between the magnitude of the capacitor voltage and 
the compensator bus as explained in the following sections. 

The matrices A, B,C,D are different for the case of ideal and 2 level inverter and are 
given in the respective sections. 

6.3 Flow chart for the process of compcnsatioii 

The flow chart for the process of compensation with an ideal current source is given in 
the following pages .The inverter logic, which is assumed ideal in the first case and in the 
second case 3phase 21evel type, is explained in their respective sections . 

The flow chart for various functions which have been caUed are also given in the 
subsequent pages. 
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mre 6.4 : Flow chart for the function update / ref 
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6.4Compensatioii with an ideal current source 
6.4.1 Introduction 


The ideal current source used for the compensation of arc furnace is assumed to have a 
constant rate of change of current whose magnitude is decided by hysteresis band and the 
maximum switching frequency allowed for the invertfer to operate. 

Ideally the current source should have an infinite band width but in our case as the 
compensator is attached to a bus where the voltages themselves are unbalanced and 
inductances are present on the either side of the compensator bus thus putting a limit on 
the maximum magnitude of the rate of change of current . Without this limit we will have 
sharp voltage spikes at the compensator bus tending to make the system unstable. 

For the purpose of our study we have fixed the following parameters for the current 
source and the system as such, 

> The rate of change of current through the current source has been fixed at IkA per 
msec. . 

> The hysteresis band has been fixed at ± 10^1 around the actual reference current. 


Thus we can say that, 

F = 

dt 

For this case the matrices B, C,D are as follows, 


^ (Rs+Rn+Rg+R) 

{Ls + ^1 + ) 


{R(2 + R) 


\ 




0 


(L^ + Lfi + Lf2 ) 
0 


6.5 
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r 


5 = 


1 


+ Zji + L^2 ) 
0 




\ 


(Z/j + + L^2 ) 

1 


^ ( {(^s ^fi)+ (Z-j + Z,j)* ^iil — {L g + Lf\) ‘ A\2) 


D = {- {(L, + Ai )• 5ii - 1} - (L, + L,i )• 5,2 ) 


6.4.2 Reference current generation and inverter control 

The system of equations is converted to discrete form and the required reference current 
is evaluated after each time step. In order to compensate for the entire reactive power and 
the alternating real power Le. (qf & p) the reference current required is given by, 


I f=I f — — V 

^ref ^af P 

where , 

Pi is the average real power taken by the load during the last half cycle . 
Thus the control flow for the current source is as follows, 


6.6 
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No change in 


< 

▼ 

Figure 6.6 : Flow chart for the current source control 

6.4.3 Arc current source current a nd their spectrum 

i 

It can be seen from the figures 6.7, 6.8 and 6.9 that we have a very significant 
improvement in the source current waveform and its harmonic content. The figures below 
show that there is a reduction in the r.m.s. cunent taken from the source because of the 
compensation of the entire reactive power. 
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The figures also indicate that the power factor has been improved to a value near unity , 
further , because ot the compensation of harmonic real power the variations in the source 
current are reduced significantly . 




Figure 6.7 : Arc and supply cuirent 



Load and Source current ( A ) supp 
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Figure 6.9b : Source current harmonic spectrum 



The reference current has been generated in such a way so as to compensate both the 
entire reactive power as well as harmonic real power. The figures 6.10 and 6.11 below 
shows the waveform of the required compensator current along with the actual 
compensator current that tracks the reference current within a band of lOA around the 
reference value . 
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Figure 6.10 : Reference and compensator current 



Reference and Compensator current { A 
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The figure below shows a blown up view of the compensator current tracking the 
reference current. 



Time ( Sec. ) 


Figure 6.11 : Compensator current tracking the reference current 

It can be seen from the figure above that infact the compensator does track the reference 
value within the specified band . The switching frequency of the compensator is variable 
with a maximum rate of 5kHz. Which is obtainable with the existing high power 
switching devices. 


6.4.5 Real and reactive power 


The figures 6.12 and 6.13 indicate the real and reactive power being consumed by the 
load and that being supplied by the source, 





71 



Time ( Sec. ) 

Figure 6.13b : Load real power 


•It is clear from the figures above that the compensated load behaves in a much better way 
than an uncompensated load because now the load draws negligible reactive power and 
less oscillating real power, thus achieving a very good compensation . 



The compensation has a boosting effect on the voltage at the critical bus. It will been 
shown below that instantaneous compensation not only improves the voltage profile but 
also reduces the harmonic content considerably. 
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Figure 6.14 : Source and critical bus voltages 

We can see from the figure 6.14 that just after the compensation begins there is a 
improvement in the voltage profile at the critical bus. Although we have some high 
frequency harmonics in the waveform but they are generally not severe on the system. 


X 10 
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Figure 6.15 : r.m.s. voltage at the critical bus 
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Comparing this Hgurc with the corresponding figure in chapter 3 we can see that the rms 

value has increased to withm ± 1% of the average rms which has increased to around 
98.5 % of the supply rrns value. 

Figure 6. 16 shows the harmonic content of the critical bus voltage, 



Figure 6.16 ; Harmonic spectrum of critical bus voltage 


Comparing this figure with the corresponding figure in chapter 3 we can see the 
tremendous improvement in the harmonic content. 

I 

The level of flicker (A)^|o ) in this case turns out to be, 

AFio=10.8 

T his value is substantially lower than that in the previous chapter and shows the dramatic 
improved of around 83% achieved with instantaneous compensation. 
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ConipcnsHtioii with h cu rrent controlled 3ph8S6 2icvcl Inverter 
6.5.1 Introduction 


In this section we will simulate a 21evel 3phase inverter and control it to act as a current 
source which will be made to track the reference current for compensation of the load. 
The circuit diagram for the inverter is given below , 



Figure 6.17 : 2Ievel 3phase inverter configuration 






75 


We have three independent single phase bridge inverters with four switehes in each. The 
switches are used in pairs to control the current through the inverter as foUows, 

Switches 1 & 2 are operated simultaneously for any phase if we are required to increase 
the current through that phase of the compensator and the switches 3 & 4 are operated 

simultaneously if we ate required to decrease the current through a particular phase of the 
compensator. 

The switching logic is dictated by the reference current being generated after each time 
step . If the compensator current is less than the lower value of the hysteresis band we 
turn on switches 1 & 2 and if the current is more than the upper value of the hysteresis 

band we turn on switches 3&4 .This logic is evaluated after each time step thus 
compensating the load instantaneously. 

The capacitor in the inverter will discharge if we do not its voltage , hence in order to 
maintain its voltage we have a control loop which generates a loss signal which is 
proportional to the voltage difference between the reference value of the capacitor 
voltage and the actual capacitor voltage . This signal is incorporated in the reference 
current generation algorithm so as to draw the requisite amount of real power jfrom the 
source and thus maintain the capacitor voltage. The control flow for the process is given 
in the next section. 


6.5.2 Reference current generation and inverter control 

The reference current required to compensate the arc furnace is generated by the 
algorithm stated in chapters. When compensating with a non ideal inverter we have to 
take into consideration the losses associated with the inverter or else the capacitor may be 
drained of charge. 
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The inverter can be modeled as a capacitor connected in series with a resistance and an 
interfecing inductance that limits the rate of rise or decay of current through the inverter 
thus checking discharging of the capacitor. 

The single phase equivalent for the inverter configuration is shown in figure 6.18, 
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Figure 6.18 : Single phase equivalent of the inverter configuration 


The differential equations governing the inverter is as follows, 


where, 


Via ~ ■ ^ca ’ 


dl 


ca 


^ib =^f 'hb+^f 


dt 

dicb 


^ic ~ * ^cc ’ 


dt 

dl 


CC 


Ict—Cf 


dt 
dVr 


dt 


+ Vp, 

6.7.1 

+ Vpb 

6.7.2 

+ Vp, 

6.7.3 


6.8 



Vja,Vii,,Vic are the inverter terminal voltages and can be either +VcOr-Vc 
depending upon the switching state at that instant . 

Ica , I cb > ^cc inverter phase currents and = ^ca ■*" ^cb ^cc • 

Vp ^ , Vpt , , Vpc are the phase voltages at the critical bus . 
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A lo=s signal representing the errer in the eapacHor voltage with respect to the reference 
voltage is generated and fed to a proportional controller. The output of the controUer, 
^loss ’ update the required reference current as follows , 


where, 



6.9 

6.10 



Figure 6.19 : Flow chart for the control of inverter 
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6 S3 Arc & source curre nt and thrir 

The figure 6.20 show the load and the source current waveforms . It can be seen from the 
figures that there is a significant improvement in the waveform of the source current as 
compared to the load current .The random fluctuations of the arc current are reduced and 
so is the rms value of the source current . 




Figure 6.20 ; Load and source current 
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The harmonic content in the current drawn from the source is reduced considerably as 
shown in the figures 6.21a and 6.21b. 

loop ^ ' r- — r 1 , , , , , 3 
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Figure 6.21a : Load current harmonic spectrum 
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Figure 6.21b : Source current harmonic spectrum 

The figures shown above indicate the effectiveness of the method of compensation using 
a voltage source inverter and the algorithm proposed by Watanabae for the instantaneous 
compensation of the entire reactive and the harmonic real power. 

The figure 6.22 below also highlight the improvement in the power factor of the arc 
furnace as seen by the source, 



6.5.4 Compensator current and reference current 

The figures 6.23 shows the required le the reference current along with the actual 
compensator current. It can be seen that the compensator current tracks the reference 
current but with a small time lag which is because of the time constant of the inverter. 
This value cannot be made very small or else we will have sharp spikes in the current 
which will be reflected in the voltage waveform also. 


< 
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6.5.5 Capacitor voltay p 


The figure 6.24 shows the waveform of the capacitor voltage of the inverter. It can be 
seen that voltage is in fact maintained through proper modification of the required 
reference cmjent taJcine into account the losses of the inverter confieuration. 
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Figure 6.24 ; Capacitor voltage of the inverter 
6.5.6 Real and reactive power 


The real and reactive being consumed by the load and that being supplied by the source is 

shown in the figures 6.25 and 6.26 given below, 

X 10 
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Figure 6.25b : Source real power 

It is clear from the above figures that the real power being supplied by the source has 
relatively lower fluctuations. The figures below show the reactive power being consumed 
by the load and that being supplied by the source, 
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Figure 6.26 : Load and source reactive power 



83 


The above figure shows the reduction in the reactive power that the source has to supply 
after the compensation of the reactive current in the load circuit. Thus it can be said that 
now the load is virtually a balanced resistive type of load to a great extent. 


6.5.7 Voltage and flicker at the critical bus 


The voltage profile at the critical bus is improved considerably after compensation . The 
figure 6.27 shows the critical bus voltage along with the supply voltage. It can be seen 
that the r.m.s. value of the bus voltage after compensation is almost equal to the supply 
voltage but some voltage spikes are introduced in the voltage waveform because of the 
rapid changes in the inverter current due to switcWng in the inverter circuit. 

The figures also indicate the reduction in the harmonic content of the voltage at the 
critical bus. 
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The figure 6.28 shows the rms value of the bus voltage. By comparing this figure with the 
corresponding figure in chapter 3 we can see the improvement in the average r.m.s. value 
after compensation. 


X 



The figure 6.29 given below shows the harmonic spectrum of the voltage at the critical 
bus. It can be seen that all the low frequency harmonics are suppressed significantly thus 
improving the quality of power at the critical bus. 
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The reduction in the harmonic content of the bus voltage suppresses the problem of 
flicker as shown below. 

The level of flicker at the critical bus , after compensation , is found to be , 


AFio=12.8 


this is an improvement of around 78% and is quite close to that obtained for the case of 
ideal current source and highlights the accuracy of the theory . for the purpose of 
instantaneous conq>ensation of arc furnace systems. 


6.6Coiiclusions 


It has been shown in this chapter that we can achieve excellent compensation by using 
the theory mentioned m chapter 5 for instantaneous compensation of different 
components of power consumed by the load. 

Thus this chapter highlights the importance of instantaneous current compensation 
scheme for effective control of the problem of flicker. It can be seen that the flicker level 
has been reduced by as much as 83% in case of compensation with an ideal current 
source and by 78% in case of compensation with a 21evel 3phase inverter. 

i 

It can also be seen that we have some spikes in the voltage waveform at the critical bus 
which are due to different rates of change of current through the inverter during different 
time intervals. These spikes can be reduced by using multilevel inverters and also by 
using different switching schemes . 
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Chapter 7 


CONCLUSIONS 


In this thesis an arc model has been constructed so as to simulate an actual arc for 
randomness, harmonic content and power utilization. Arc furnace installation has been 
simulated and its behavior highlighted without any compensating aids thus bringing out 
the necessity of compensation. Flicker has been defined and evaluated at the critical 
bus.It has been shown that the arc current has a large harmonic content and that the arc 
furnace draws randomly fluctuating real and reactive powers from the supply, leading to 
voltage fluctuations at the critical bus .A high flicker level brings out the importance of 
flicker management. 

Next we simulated an arc furnace system with a shunt connected S.V.C. It has been 
shown that although we are able to bring down the average reactive power required by 
the furnace but still a S.V.C as such is incapable of reducing the harmonics to any 
considerable extent. It has also been shown that a S.V.C has a limited affect as far as 
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flicker IS concerned basically because it compensates for the fundamental reactive current 
only. Thus, it has been brought out that an arc fiimace does require a fast and rapidly 
variable compensation so as to put a check on the problem of flicker. 

Knowing the need of fast compensation, various instantaneous compensation schemes 
were considered and a new scheme of compensation based on the theory proposed by Lai 
[10] and that proposed by Watanabae [13] has been considered. 

Next, the same arc furnace system has been simulated, but with the new compensation 
scheme using an ideal current source. It has been shown that the harmonic content in the 
supply current and voltage at the critical bus were reduced to a great extent. The reactive 
power required from the supply also went down to near zero and also there was a 
reduction in the fluctuations of real power. Flicker was reduced substantially thus 
showing a excellent improvement in the behavior of the system. Then we realized the 
ideal current source with a 21evel 3phase inverter using current feedback. It has been 
shown that the system still behaves considerably better than those that have been 
compensated using S.V.C’s. 

Thus it was brought out that instantaneous current compensation method described in this 
study has excellent performance as compared to conventional methods. Although the 
21evel 3phase inverter introduces harmonics at switching frequency but still it reduces 
flicker to substantially lower levels. 

Scope for further woik 

The following the areas which still require some attention, 

> The inverter used for compensation needs to be upgraded in the sense that the 
switching frequency should be made more or less constant so that we have lesser 
spikes in the current waveform which are reflected in the voltage at the critical bus. 
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^ Practical implementation of the theory could be done so as to verify the theoretical 
results obtained in this study. 

^ Investigations should be done so as to see the affect of replacing the 21evel inverter 
with higher level inverters. 

> Modifications in the reference current generator scheme should be investigated so as 
to sense the fluctuations in the real power drawn by the arc furnace more accurately. 
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APPENDIX 


A.1 Program for compensation of are furniice with S V.r. 


Clear all 


global 

k 

Td 

t 





global 

w 

vs 

Va 

Vb 

Vc 



global 

Rs 

Lsl 

Ls2 

Ls 

Rt 

Lt 


global 

R1 

R2 

R3 

R4 

R5 

Ra 

Rb Rc 

global 

Lf 

Cf 

L 

Ilo 

Ico 

Ka 

Kb Kc 

global 

ALFa 

ALFb 

ALFc 

ja 

jb 

jc 

ThetaA ThetaB ThetaC 

global 

Aoa 

Boa 

Aob 

Bob 

Aoc 

Boc 


global 

Ala 

Bla 

Alb 

Bib 

Ale 

Blc 


global 

A2a 

B2a 

A2b 

B2b 

A2c 

B2c 


global 

Xa 

Xb 

Xc 

Q1 

Qs 

Vp 


global 

Via 

V2a 

Vlb 

V2b 

Vic 

V2c 



Tdl-1/36000; 

w - 100*pi ;V3-33000*sqrt{2/3) ; 

Rs- 0.254 ; Lsl-O. 00006 ; 

Rt-0.06 ; LS2-0. 00238 

Ls « Lai + L32 ; Lt - 0.03 ; 

Rl-3 ; R2-8 ; R3-10 ; R4-15 ; R5 - 20 ; 

C£ « 0.0001 ; L - 0.05 ;Lf - l/(9*w*w^Cf) ; 

Ilo - Vs/(w’^L) ; ICO Vs*w*Cf/(w*w*Lf*Cf-l) ; 
Tc * 0.1 ; Tf - 0.5 ; 

k « 1 ; 
tU) - 0 ; 

Va(k) - Vs^sintw’^t (k) ) ; 

Vb(k) * vs*sin(w*t (k) - 2*pi/3) / 

Vc{k) « V»*sin{w*t (k) + 2*pi/3) ; 


Xa(:,l) - [ 

0 

; 

0 

; 0 / 0 3 

/ 

Xb(:,l> - [ 

0 

; 

0 

; 0 ; 0 ] 

; 

Xc ( : , 1 ) - [ 

0 

; 

0 

; 0 / 0 3 

; 

Vp(:,l) - ( 

0 

; 

0 

; 0 3 ; 


Ql(:,l) - 1 

0 

; 

0 

; 0 3 ; 


Qa(;,l) - C 

0 

; 

0 

; 0 3 ; 


ja - 0 ; jb 


0 

# 

jc - 0 ; 


Ka - 0 ; Kb 

M 

0 

; 

Kc - 0 ; 



y j | y ^ 

ThetaA - 4000 ; ThetaB - 4000 ; ThetaC - 4000 ; 

Ca - 0 ; Cb - 0 ; Cc - 0 ; 

while { t{k) < Tc ) 

k - k + 1 ; 

Matrix ; 

Xa(:rk) * uncomp ( Aoa,Boa,Xa(:/ (k-1) ) /Va(k*-1) ) ; ja * ja + 1 r 

Vp(l^k) ■ uncompBusV( Xa(2, (k-1) ) #Ra,Va(k“l) ) ; 

Xb(S|>k} *• uncomp ( Aob^Bob/Xb ( i p (k~l) ) fVb(k""l) ) <- jb jb + i , 

Vp(2,k) * uncompBusV( Xb(2, (k-1) ) #Rl>,Vb(k-l) ) ; 

Xc(:,k) *- uncomp ( Aoc,Boc,Xc(:f (k-1) ) ,Vc(k-l) ) ; jc - jc + 1 , 

Vpd^k) * uncompBusV( Xc (2, (k~l) ) /Rc,V'c(k-l) ) ; 


updateTVQ ; 

end 


ThyA ** 
while 


0 ; ThyB « 0 ; ThyC - 0 ; , v x 

( ( t(k) >- TC )fi{ t(k) < (TC + 2*pi/w) ) 
k - k + 1 ; 


) 


if 


Matrix ; 

( (peakvoltage ( Vpd/S) 


)) I (Ca 


1 ) ) 
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Ca - 1 ; 

Xa(:,k) - cpcomp( Ala,Bla,Xa( (k-i) ) va(t n » 

alaait 7cl'"l ^^«'-2a,Ka,Xa(:/,k-l„,k’(k-l, 


ja + 1 


el3© 

end 

if 


Xa(:,k) - uncorap( Aoa,Boa,Xa( i , (k-l) ) ,va(k-ii ^ x i 

Vp(l,k) - uncompBusV( xa(2, (k-l)),Ra,Va k-1 ja+ 1 ; 

error (' unexpected situation') ; > '• 


( (peakvoltage( vp(2,:) ) ) | (Cb 

Cb - 1 ; 


1 ) ) 


Xb(:,k) - cpcomp( Alb,Blb,Xb ( i, (k-1) ) , vb(k-l) ) • ib - ih + i ■ 


else 

end 

it 


Aob,Bob,Xb{!, (k-l) ),vb(k-l) ) ;jb = jb + 1 ; 
Vp(2,k) - uncompBu3V( Xb(2, (k-1) ) ,Rb,Vb(k-l) ) ; 

error (* unexpected situation’) ; 


end 


( (peakvoltage( Vp(3,j) ) ) | (Cc — l ) ) 

Cc - 1 ; 

Xc(:,k) - cpcomp( A1 c,B1c,Xc( : , (k-1) ) ,Vc(k-l) ) ;jc - jc + 1 ; 
Vp(3,k) - oompBu8V( vic,V2c,Rc,Xc(:,(k-l)),Vc(k-l),ThyC ) ; 

elseif { CC 0 ) r f 

Xc(i,k) * uncompl Aoc,Boc,Xc ( j , (k-1) ) , Vc(k-l) ) ;jc - jc + 1 ; 
Vp(3,k) - uncompBu8V( xc(2, (k-l)),Rc,Vc(k-l) ) ; 

els© error ( ’unexpected situation') ; 
end 
updateTVQ / 


ThyA • 0 / ThyB * 0 ; ThyC • 0 ; 

AbortA - 0 ; AbortB « 0 ; AbortC - 0 ; 
while ( t(k) < Tf ) 
k - k + 1 ; 
if { k < ThetaA ) 

if ( rerocroflsingi Vp(l,i) ) ) 

ThyA « 0 t 
AbortA « 1 I 
el«e ThyA *■ 0 } 
end 

elaeif ( k •••* ThetaA ) 

if ( AbortA — 0 ) 

ThyA “ 1 ; 

elaeif ( AbortA — 1 ) 

ThyA - 0 f 

else error {’check firing’) ; 

end 

elaeif ( k > ThetaA ) 

if ( AbortA — 0 ) 

if {( zerocrossing ( Xa(4,:) ) — 0 )&( ThyA — 1 )) 

ThyA * 1 ; 

elseif (( zerocrossing ( Xa(4,:)) — 1 )l( Xa(4,(k-1)) 0 ) ) 

ThyA » 0 ; ■ 

Xa(4, (k-D) - 0 / 

end 

elaeif ( AbortA »• 1 ) i 

ThyA - 0 ; 

else error (’check after firing instant of ThyA’) ? 
end 

else error (* unexpected logic*) ; 
end 

if ( k < ThetaB ) 

if ( zerocrosfling ( Vp(2,:) ) ) 

ThyB - 0 ; 

AbortB - 1 ; 
else ThyB - 0 t 
end 

elseif ( k ThetaB ) 

if ( AbortB — 0 ) 

ThyB • 1 ; 

elseif ( AbortB 1 ) 

ThyB - 0 / 
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else error ( 'check firing') ; 

nnd 

( k > ThetaB ) 

if ( Aborts -« 0 ) 

if (( zerocrossing{ Xb(4,:) ) — O )&( ThyB — ' 
ThyB - 1 ; 

elseif (( zerocrossing ( Xb(4,:)) — 1 )|( Xb(4, 
ThyB - 0 ; 

Xb(4, (k-D) » 0 ; 

end 

elseif ( Aborts 1 ) 

ThyB - 0 ; 

else error (* check after firing instant of ThyA') ; 
end 

else error ( ’unexpected logic*) ; 
end 

if ( k < ThetaC ) 

if ( zerocrossing { Vp{3,:) ) ) 

ThyC * 0 ; 

AbortC * 1 ; 
else ThyC - 0 ; 
end 

elseif ( k ThetaC ) 

if ( AbortC 0 ) 

ThyC - 1 ; 

elseif ( AbortC •— 1 ) 

ThyC - 0 ; 

else error (* check firing’) ; 
end 

elseif ( k > ThetaC ) 

if ( AbortC — 0 ) 

if (( zerocrossing ( Xc(4#.:) ) 0 )&( ThyC ** 

ThyC - 1 ; 

elseif (( zerocrossing ( Xc(4,:)) — 1 )|( Xc(4, 
ThyC - 0 ; 

Xc(4, (k-D) - 0 ; 

end 

elseif ( AbortC -- 1 ) 

ThyC » 0 ; 

else error (’check after firing instant of ThyA’) ; 
end 

else error ( ’unexpected logic') ; 
end 


Matrix ; 

X«(j,k) - flcomp(Ala,Bla,A2a,B2a,Xa(:, (k-D), Va(k-l), ThyA) ;ja 
Vp(l,k) - compBusV( Via, V2a,Ra,Xa(:, (k-D), Va(k-l), ThyA ) ; 

Xb(;,k) - flcoinp(Alb, Bib, A2b,B2b,Xb(:, (k-D), Vb(k-D, ThyB) ;jb 
Vp(2,k) - compBu3V( Vlb,V2b,Rb,Xb(:, (k-D ),Vb (k-D, ThyB ) ; 

Xc(:,k) - flcomp(Alc,Blc,A2c,B2c,Xc(:, (k-D), Vc(k-D, ThyC) ;jc 
Vp(3,k) - coinpBusV( Vic, V2c,Rc,Xc(:, (k-D ),Vc(k-D, ThyC ) ; 

' I 

updateTVQ ; 


if ( ( ThyA — 0 ) & (peakvoltage { Vp(l,;) ) — 1) ) 
Iqa - updatelQ( 01(1, i) ) ; 

ThetaA « updateALFA(Iqa) ; 

AbortA - 0 ; 

end 

if ( ( ThyB — 0 )& (peakvoltage ( Vp(2,;) ) «» 1) ) 

Iqb - updateIQ( Ql(2, :) ) ; 

ThetaB « updateALFA(Iqb) ; 

Aborts » 0 ; 

end 

if ( ( ThyC — 0 )& (peakvoltage ( Vp(3,:) ) — 1) ) 
Iqc - updateIQ( Ql(3^:) ) ; 

ThetaC * updateALFA(Iqc) ; 


I )) 

(k-D) — 0 )) 


1 )) 

(k-D) — 0 )) 


ja + 1 ; 
jb + 1 ; 

3C + 1 ; 


end 


AbortC »« 0 ; 


end 


%%%%%%%%%%»%%%%%%%%%%% 


function Vp - updateBusV( ?! , R , X , 0 ) 
global Lt Rt 

Vp - Z(l)*Lt*X(l) + (Lt*Z(2)+R+Rt)*X(2)+Lt*Z(3)*X(3)+Lt*Z(4)*X(4)+Lt*Z(5)*U ; 

function RTN - updateALFA(Iq) 
global k Ilo Ico 

alf(l) * pi/6 ; itemp(l)-0 ; Icount (1) »*0 / 
fire - pi/2 ; 
for ( p-2:180 ) 

alf(p) - alf(p~l) + pi/540 ; 

itemp(p) - 2*Ilo*( pi/2 - alf(p) - (sin {2-^alf (p) ) ) /2 ) /pi ; 
Icount(p) - Ico - iteinp(p) ; 
if ( Icount(2) > Iq ) 
fire - pi/6 ; 

elseif (( Icount(p) > Iq )&( Icount(p-l) < Iq ) ) 
fire » alf(p) ; 

end 

end 

RTN - round (fire^360/pi 4- k) ; 

%%%%%%%%%«%%%%%%%%%%%%%%«%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
function Vol - uncompBusV( I , R , U ) 

global Rs Ls Rt Lt Lsl Ls2 

Vp - (Ls^(R+Rt)-Rs*Lt)*I/(Ls+Lt) + Lt*U/(Lt+Ls) ; 

Vol - Vp ; 


111111%% III %%%!%% !%%%%%! %%%%%%%%%!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%% 
function state • uncomp (F,G/X,U) 

X - F*X 4- G*U ; 

State - X ; 

III I II III III! I III! Ill III 11 mill iiiiiiiiiiiiiiiiii 
function peak - peakvoltage (V) 

global k 

peak ■* 0 ; 

if ( ( aba(V(k-l)) < aba (VOc-Z) ) ) i ( ab8(V(Ic-2)> > ab8(V()c-3)) ) ) 
if ( V(k-l) > 0 ) , 

peak - 1 ; 

else if { V(k-l) < 0 ) 
peak ** -1 ; 

end 

end 

peak * abs (peak) ; 

else 

peak ■“ 0 ; 
end , 

function Matrix 

global ja Ka Ra jb Kb Rb jc Kc Rc 
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if ( ja >« Ka ) 
ja « 0 ; 

Tamp - updateR ; 

Ra -• Temp(l) ; 

Ka - Temp (2) ; 
updatematA ; 

end 

if ( jb >» Kb ) 
jb « 0 ; 

Temp * updateR ; 

Rb « Temp(l) ; Kb - Temp(2) ; 
updatematB ; 

end 

if ( jc >- Kc ) 
jc - 0 ; 

Temp - UpdateR ; 

Rc - Temp(l) ; Kc - Temp (2) ; 
updatematC ; 

end 

function state - flcomp(Fl,Gl, F2,G2,X,U,Thy) 

if ( Thy — 0 ) 

X * Fl’^X + G1*U / 

©Iseif ( Thy -* 1 ) 

X - F2*X 4- G2^U ; 

end 

state X ; 

function state « cpcomp (F, G,X,U) 

X - F^'X + ; 

state - X ; 

function Vp « compBusV( Z1 , 22 , R , X , U ,Thy ) 
global Lt Rt 
if { Thy — 0 ) 

Vp - Zl{l)‘^Iit'^Xa) + (Lt*Zl(2)-fR4<Rt)^X(2)+Lt*Zl(3)*X(3)4'Lt*Zl(4)*X(4)+Lt*Zl(5)*U ; 
elseif ( Thy — 1 ) 

Vp *• Z2(l)*Lt*X(l) + (Lt*Z2(2)+R4*Rt)*X(2)+Lt*Z2(3)*X(3)+Lt*Z2(4)*X(4)+Lt*Z2(5)*U ; 
end 

1 111% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%% 
function RTNI - updatelQ(Q) 

global k Vs 

S » 0 ; I 

if ( k > 360 ) 

for m * 1:360 i 

S - S + Q<k+l-m) ; 
end I 

end 

Qavg *- S/360 ; 

RTNI « Qavg/Vs ; 

111111%%% III 

function cross - zerocrossing ( Y ) 


cross ** 0 ; 

if ( ( Y(k-2) < 0 )&( Y(k-l) > 0 ) ) 
cross « 1 ; 

elseif ( ( Y(k-2) > 0 )&( Y(k-l) < 0 ) ) 
cross *" “1 ; 


global k 
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else cross » 0 ; 

end / 

cross - abs (cross) ; 

function updateTVQ 

global k t Td w Vs Va Vb VC Xa Xb Xc Qla Qlb Qlc Qsa Qsb Qsc Vp Q1 Qs 
global Vpa Vpb Vpc 

t(k) « t(k"l) + Td ; 

Va(k) - Vs*sin(w*t(k)) ; 

Vb(k) - Vs*sin(w*t(k) - 2’^pi/3) / 

Vc{k) « V8*8in{w*t(k) + 2*pi/3) ; 

Ql(l,k) « -( ( Vp(2,k) * Xc(2,k) ) - ( Vp(3,k) * Xb(2,k) ) ) ; 

Ql(2,k) - -( ( Vp(3,k) * Xa(2,k) ) - ( Vp(l,k) * Xc(2,k) ) ) ; 

Ql(3,k) - -( ( Vp(l,k) * Xb(2,k) ) - ( Vp(2,k) * Xa(2,k) ) ) ; 

Qs{l,k)*-( (Vb(l,k)’^(Xc(2,k)+Xc(3,k)+Xc(4,k)))-(Vc(l,k)*{Xb(2,k)+Xb(3,k)+Xb(4,k)))) ; 
Q8(2,k) — ( (Vc(l,k)*(Xa(2,k)'fXa(3,k)4-Xa(4,k)))-(Va(l,k>’^(Xc(2,k)fXc(3,k)fXc(4,k))))/ 
Qa(3,k)*— •{(Va(l,k)^(Xb(2,k)+Xb(3,k)fXb(4,k)))-(Vb(l,k)^(Xa(2,k)fXa(3,k)+Xa(4,k)))); 

%% 9 ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%«%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%% 
function X updateR 

global R1 R2 R3 R4 R5 

m - rand* 10 ; 

if ( (m>-0)&(m<2) ) 

R - R1 ; K “ round (rand* 100) ; 
eiseif ( (in>*2 ) & (ni<4 ) ) 

R R2 ; K - round (rand* 3 00) ; 
eiseif ((m>«4)&(m<6)) 

R » R3 ; K » round (rand* 2 00) ; 
eiseif ( (ra>*6) & (m<8) ) 

R - R4 ; K - round (rand* 100) / 
eiseif ( (m>-8) & (m<10) ) 

R - R5 ; K * round (rand* 50) f 

end 

X - [ R K 3 ; 

function updatematA 
global Td 

global Rs Ls Rt Lt Ra 
global Lf Cf L 

global Aoa Boa Ala Bla A2a B2a 
global Via V2a 
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al - 

( 0 


0 


1/Cf 


0 



Ml 


Mi! 


M3 ' 


0 



(Lt*Ml - 1)/Lf (Lt*M2 

+ (RtRt); 

l/Lf 

Lt*M3/Lf 


0 



0 


0 


0 


0 

] ; 

bl « 

[ 0 


M5 

t 

Lt*M5/Lf 

/ 

0 

] ; 

a2 - 

[ 0 


0 


1/Cf 


0 



N1 


N2 


N3 


N4 



(Lt*Nl - 1)/Lf (Lt*N2 

+ (R+Rt) 

)/Lf 

Lt*N3/Lf 


Lt*N4/Lf 



Lt*Nl/L (Lt* 

N2 

+ (R+Rt))/L 

Lt*N3/L 


Lt’*^N4/L 

] ; 

b2 - 

C 0 


N5 

f 

Lt*N5/Lf 


Lt*N5/L 

] ; 

[Aoa, 

Boa] - c2d{ao,bo,Td) 

; 







[Ala, 

Bla] - c2d(al,bl,Td) 

} 







[A2a, 

B2a] - c2d(a2,b2,Td) 

t 







Via - 

■ [ Ml ; M2 ; M3 ; M4 

; 

M5 ] ; 






V2a - 

• [ N1 ; N2 ; N3 ; 114 

; 

N5 ] ; 







%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
function updatematB 


global Td 

global Rs Ls Rt Lt Rb 
global Lf Cf L 

global Aob Bob Alb Bib A2b B2b 
global Vlb V2b 


R - Rb ; 

DM « Ls*Lt*(l/Lf + 1/Ls + 1/Lt) ; 

DN - Ls*Lt*(l/Lf + 1/Ls + 1/Lt + 1/L) ; 

Ml « Ls/<Lf’^DM) ; N1 

M2 « -(Rs+(LS/Lf + 1)*(R+Rt) )/DM ; N2 

M3 - -Rs/DM ; N3 

M4 « 0 ; N4 

M5 - 1/DM ; N5 


Ls/(Lf*DN) 

-(Rs+(Ls/Lf + LS/L -f 1)*(R+Rt) )/DN 
-Rs/DN 
-Rs/DN 
1/DN 


ao 


bo 

al 


0 

0 

0 

0 

0 

0 

Ml 


-(R+Rt+Rs)/(Lt4-Ls) 

0 

0 

l/(Lt+Ls) 

0 

M2 


(Lt*Ml - 1)/Lf (Lt^M2 + (RtRt))/Lf 


bl 

a2 


0 

0 

0 

N1 


0 

M5 

0 

N2 


(Lt^Nl - 1)/Lf (Lt*N2 + (R+Rt))/Lf 


0 

0 

0 

0 

0 

1/Cf 

M3 

Lt*M3/Lf 

0 

Lt*M5/Lf 

1/Cf 

N3 

Lt^N3/Lf 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

N4 

Lt*N4/Lf 


] f 
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Lt*Nl/L (Lt*N2 + (R+Rt))/L Lt*N3/L Lt*N4/L ] ; 

t 0 '■ N5 ; Lt*N5/Lf ; Lt*N5/L ] ; 

[Aob,Bob] - c2d(ao,bo,Tcl) ; 

[Alb.Blbl - c2d(al,bl,Td) ; 

[A2b,B2b] - c2d(a2,b2,Td) ; 

Vlb - [ Ml ; M2 ; M3 ; M4 ; M5 ] ; 

V2b - I N1 ; N2 ; N3 ; N4 ; N5 ] ; 

%%%%%%%%%*%%%»%%%%%%%%%%%»%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function updatematc 

global Td 

global Ra La Rt Lt Rc 
global Lf Cf L 

global Aoc Boc Ale Blc A2c B2c 
global Vic V2c 

R - Rc ; 

DM - L3*Lt*(l/Lf + 1/Ls t 1/Lt) ; 

ON - Ls*Lt*(l/Lf + 1/Ls + 1/Lt + 1/L) ; 


Ml » 

Ls/{Lf*DM) 

; Nl 

* Ls/(Lf*DN) 


M2 - 

-(R8+{Ls/Lf + 1)* 

(R+Rt))/DM ; N2 

- -(Rs+(Ls/Lf 

+ Ls/L + 1)*(R+Rt))/DN 

M3 - 

-RS/DM 

; N3 

=« -Rs/DN 


M4 • 

0 

; N4 

* -Ra/DN 


M5 - 

1/DM 

; N5 

- 1/DN 


ao * 

[ 0 

0 

0 

0 


0 

<-{R+Rt+Rs)/(Lt+Ls) 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 ] ; 

bo - 

C 0 

l/(Lt+Ls) ; 

0 

0 ] ; 

al - 

C 0 

0 

1/Cf 

0 


Ml 

M2 

M3 

0 


(Lt*Ml - 1)/Lf 

(Lt*^M2 + (R+Rt))/Lf 

Lt*M3/Lf 

0 


0 

0 

0 

0 ] ; 

bl - 

[ 0 ; 

M5 ; 

Lt*H5/ht 

; 0 ] ; 

a2 - 

C 0 

0 

1/Cf 

0 


Nl 

N2 

N3 

N4 


(Lt*Nl “ 1)/Lf 

(Lt*N2 + (R+Rt))/Lf 

Lt*N3/Lf 

Lt*N4/Lf 


Lf*'Nl/L 

(Lt’^N2 + (R+Rt))/L 

Lt^N3/L 

‘ Lt*N4/L ] ; 

b2 - 

[ 0 

N5 ; 

Lf*^N5/Lf 

; Lt*N5/L ] / 


[AOC, Boc] - c2d(ao,bo,Td) ; 

[Ale, Blc] - c2d(al,bl,Td) ; 
[A2c,B2c] - c2d{a2,b2,Td) ; 

Vic - [ Ml ; M2 ; M3 ; M4 ; M5 ] ; 
V2c - [ Nl ; N2 ; N3 ; N4 ; N5 ] ; 
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A.2 Program for compensation of arc furnace with ideal current source 


clear all 


global R1 R2 R3 R4 R5 Ro 

global Rs Ls Rtl Ltl Rt2 Lt2 Ra Rb Rc 

global ja Ka jb Kb jc Kc 

global Fla Gla Fib Gib Flc Glc 

global Ala A2a Bla B2a Alb Bib A2b B2b Ale Blc A2c B2c 

global Xa Xb Xc 

global II Fc Iq h pi Th PI 

global k t Td Tc w Vs Vo Ua Ub Uc Vp slope 

Tc - 0.3 ; Tf - 1.5 ; Td « 2/(360^60) ; 

f - 60 ; w « 2'^pi*f ; Vo « 33000*sqrt (2/3) ; 

slope 1000000 ; h - 10 ; Th - roundd/ (2*Td*f ) ) ; 

R1 - 3 ; R2 - 8 ; R3 - 10 / R4 - 15 ; R5 - 20 ; 

RS - 0.254 ; Ls - 0.00005 ; 

Rtl* 0.06 ; Ltl- 0.001 

Rt2« 0.001 ; Lt2- 0.03 


k -1 ; t(k) - 0 ; 

Vs(l/k) - Vo*sin(w^t (k) ) ; 

Vs(2^k) - Vo*8in(w*t (k) - 2*pi/3) ; 
V3(3,k) - Vorsin (w^t(k) + 2*pi/3) / 
Xa(:,k) - [ 0 ; 0 ] ; 

Xb(:,k) - [ 0 ; 0 ] ; 

Xc(:/k) - C 0 ; 0 ] ; 

Il{:,k) - C 0 / 0 ; 0 ] ; 

Is(:,k)-[0;0/0]; 

Iq( ; , k) - C 0 ; 0 / 0 ] / 

pl(l) - 0 ; Pl(l) - 0 ;ps(l) - 0 ; 

Fc(:,)c) - [ 0 ; 0 ; 0 ] ; 

Ql(:,lc)-C0;0;0]; 
Q3(:,k)-[0;0;0]; 
ja - 0 ; jb » 0 ; jc - 0 ; 

Ka - 0 ; Kb - 0 ; Kc “ 0 ; 

Ra-0;Rb"0;Rc“0; 

Ua(:,k) - ( vs(l,k) ; Fc(l,k) ] ; 

Ob(i,k) - [ V8{2,k) ; Fc(2,k) ] ; 

Uc{:,k) - [ V3(3,k) ; Fc(3,k) ] ; 


while ( t(k) < Tf ) 


k - k + 1 ; 
t(k) - t(k-l) + 
solve_phaseA ; 
solve_pha8eB ; 
solve phaseC ; 
Il(l,k) 

11(2, k) 

11(3, k) 

Is(l,k) 

Is(2,k) 

Is(3,k) 

Vp (1, k) 
vp(2,k) 

Vp(3,k) 

V3(l,k) 


Td 


V8(2,k) 
Vs(3,k) 
Ql(:#k) 
Qs ( : , k) 
lq(:,k) 

P8(k) 


Xa(l,k) + Xa(2,k) 

Xb(l,k) + Xb(2,k) 

Xc(l,k) + Xc(2,k) 

Xa(l,k) ; 

Xb(l,k) ; 

Xc(l,k) ; 

Ala*Xa(l, (k-1) ) + 

Alb*Xb(l, (k-D) + 

Alc*Xc(l, (k-1) ) + 
Vo*sin(w*t(k) ) ; 

Vo*8in(w*t(k) - 2*pi/3) ; 
Vo*sin(w*t(k) + 2*pi/3) ; 
updateQ( Vp(:,k) , Il(!,k) 
updateQ( Vs(:,k) , Is(!,k) 
updatelqC Vp(i,k) , II (’/k) 


A2a*Xa(2, (k-D) 
A2b*Xb(2, (k-D) 
A2c*Xc(2, (k-D ) 


) ; 


) ; 


- Vs (1, k) *18 (1, k) + Vs(2,k)*l8(2,k) 


+ Bla*V3(l, (k-D) + B2a*Fc(l, (k-D) 
+ Blb*vs(2, (k-D) + B2b*Fc(2, (k-D) 
+ Blc*Vs(3, (k-D) + B2c*Fc(3, (k-D) 


+ V8(3,k)*l8(3,k) ; 


if ( t(k) >= 1/60 ) 

S » 0 ; 
for j ■* l:Th 
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a - a + pti(k^j+i.) / 

Rild 

Ps(k) - S/Th ; 

•nd 

if ( t (k) > To ) 

if ( Xa(2,k) <- ( lq(l,k) - h ) ) 

Fc(l,k) “ slope ; 

elseif ( Xa{2,k) >- ( Iq(l,k) + h ) ) 

Fc(l,k) - -slope ; 
else Fc(l,k) - Fed, (k-1)) ; 
end 

if ( Xb(2,k) <- ( Iq(2,k) + h )- ) 

Fc(2,k) - Slope ; 

elseif ( Xb(2,k) >- ( Iq(2,k) + h ) ) 

Fc(2,k) “ -slope ; 
else Fc(2,k) - Fc(2, (k-1)) ; 
end 

if ( Xc(2,k) <- ( Iq(3,k) + h ) ) 

Fc(3,k) - slope ; 

elseif ( Xc(2,k) >- ( Iq(3,k) + h ) ) 

Fc(3,k) " -slope ! 
else Fc(3,k) - Fo(3,(k-l)) ; 
end 

els© Fc(:,k) - [ 0 ; 0 ; 0 ] ; 

end ' 

Ua(:,k) - [ Vs(l,k) ; Fc(l,k) ] ; 

Ob(:,lc) - [ V3(2,k) ; Fc(2,k) ] t 

Uc(;,k) - [ V8(3,k) ; Fc(3,k) ] ; 

end 

»»%»%»%%%%%%%%%%%%%%»%%%%%%« %%%%%%%%%%%%%%%%%%»%%%»%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
function aolvejphaseC 

global k Td R1 R2 R3 R4 Ro 

global Ra Ls Rtl Ltl Rt2 Lt2 Rc 

global jc Kc 

global Flc Glc 

global Xc Uc 

global Ale Blc A2c B2c 

if ( jc >- Kc ) 

jc “ 0 ;Teinp ■“ updateR ; Rc - Temp(l) ; Kc “ Tenip(2) ; 

Ale - -(RS + Rtl) + (Ls + Ltl)*(Rs + Rtl + Rt2 + Rc)/(Ls + Ltl + Lt2) ; 

A2c - (Ls + Ltl)*(Rt2 + RC)/(LS + Ltl + Lt2) ; 

Blc - - (Ls + Ltl)/(L3 + Ltl + Lt2) + 1 ; 

B2c - (Ls + Ltl)*Lt2/(Ls + Ltl + Lt2) ; 

All - -(Rs + Rtl + Rt2 + Rc)/(Ls + Ltl + Lt2) / 

A12 -• -(Rt2 + Rc)/(Ls + Ltl + Lt2) ; 

A21 - 0 ; 

A22 - 0 ; 

Bll - l/(Ls + Ltl + Lt2) ; I 

B12 - -Lt2/(LS + Ltl + Lt2) } 

B21 - 0 ; : 

B22 - 1 ; 

AC - t All A12 

A21 A22 ) ; 

Be - [ Bll B12 

B21 B22 ] ; 

[Flc Glc] - c2cl(Ac,Bc,Td) ; 

end 

Xc(:,k) - Flc*Xc(!, (k-1) ) + Glc*0c ( : , (k-1) ) ; 
jc - jc + 1 ; 
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function solvejphaseB 

global k Td R1 R2 R3 R4 Ro 

global Rs Ls Rtl Ltl Rt2 Lt2 Rb 

global jb Kb 

global Fib Gib 

global Xb Ub 

global Alb A2b Bib B2b 

if ( jb >- Kb ) 

jb * 0 ; Temp « updateR ; Rb = Temp(l) ; Kb * Temp (2) ; 

Alb « -{Rs + Rtl) + (Ls + Ltl)*{Rs + Rtl + Rt2 + Rb) / (Ls + Ltl + Lt2) ; 

A2b » (Ls + Ltl)*(Rt2 + Rb)/(L 3 + Ltl + Lt2) ; 

Bib - -(Ls + Ltl)/(LS + Ltl + Lt2) + 1 ; 

B2b - (Ls + Ltl)’^Lt2/(Ls + Ltl + Lt2) ; 

All * -(Rs + Rtl + Rt2 + Rb)/(Ls + Ltl + Lt2) ; 

A12 « -(Rt2 + Rb)/(Ls + Ltl + Lt2) ; 

A21 « 0 ; 

A22 - 0 ; 


Bll - l/{Ls + Ltl + Lt2) ; 

B12 - -Lt2/(Ls + Ltl + Lt2) ; 


B21 - 0 ; 
B22 - 1 ; 


Ab « C All 

A12 

A21 

A22 ] ; 

Bb - [ Bll 

B12 

B21 

B22 ] ; 

[Fib Gib] » 

c2d(Ab,Bb,Td) ; 


«nd 

Xb(i,k) - Flb*Xb(:, (k-D) + Glb*Ub( ! , (k-1) ) ; 
jb - jb + 1 ; 


1 %%%%%%% %%%%%%%%%!%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%• 
function solve_phaseA 

global k Td R1 R2 R3 R4 Ro 

global RS Ls Rtl Ltl Rt2 Lt2 Ra 

global ja Ka 

global Fla Gla 

global Xa Ua 

global Ala A2a Bla B2a 

if ( ja >- Ka ) 

ja - 0 ; Temp - updateR ; Ra - Temp(l) ; Ka - Temp (2) ; 

Ala - -(Rs Rtl) + (Ls + Ltl)'^(Rs + Rtl + Rt2 + Ra)/(Ls ,+ Ltl + Lt2) ; 

A2a » (Ls ’¥ Ltl)*(Rt2 + Ra)/(Ls + Ltl + Lt2) ; 

Bla - -(LS + Ltl)/(Ls + Ltl + Lt2) + 1 ; 

B2a « (Ls + Ltl)*Lt2/(Ls + Ltl + Lt2) ; 

All » -(Rs + Rtl + Rt2 + Ra)/(Ls + Ltl + Lt2) ; 

A12 - -(Rt2 f Ra)/(Ls + Ltl + Lt2) ; 

A21 -* 0 ; 

A22 « 0 ; 

Bll - l/(Ls 4* Ltl t Lt2) ; 

B12 * -Lt2/(LS + Ltl + Lt2) ; 

B21 - 0 ; 

B22 - 1 / 

Aa * (All A12- 

A21 A22 ] ; 
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Ba * [ Bll B12 

B21 B22 ] ; 

[Fla Gla] c2d{Aa , Ba , Td) ; 


end 

Xa(:,k) - Fla*Xa(:, (k-1) ) + Gla*Ua ( : , (k-1) ) 
ja - ja + 1 ; 

function X - updateR 
global R1 R2 R3 R4 R5 Ro 

m « rand’^lO ; ^ 

if ( (m>-0) & (ni<2) ) 

R - R1 ; 

elseif ( (m>»»2)r& (m<4) ) 

R « R2 ; 

elseif ( (in>-4) & (m<6) ) 

R - R3 ; 

elseif { (m>»«6) & (m<8) ) 

R * R4 ; 

elseif ( (m>-8) & (m<10) ) 

R » R5 ; 

end 

if ( R — R1 ) 

K - round (rand* 100) / 
elseif ( R — R2 ) 

K -« round (rand*3 00) ; 
elseif ( R R3 ) 

K - round (rand*2 00) ; 
elseif ( R — ' R4 ) 

K - round (rand* 100) ; 
elseif ( R R5 ) 

K * round (rand*50) ; 

end 

X - [ R K ] ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
function Q * updateQ.( V , I ) 

ql - V(2)*I(3) - V(3)*I(2) ; 
q2 - V(3)*I(1) - V(1)*I(3) ; 
q3 - V(1)*I(2) - V(2)*I(1) ; 

Q - [ -ql / -q2 ; -q3 ] ; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
function P - update? ( V ^ I ) 

P » V(1)*I(1) + V(2)*I(2) + V(3)*I(3) ; 

%l%l %%%%%%% 1%%%%%%%%%%%%%%%%%%%%%% 
function Current “ updatelq( V , I ) 

global t k pi Th PI 

NormV - V(1)*V(1) + V(2)*V(2) + V(3)*V(3) ; 
pl(k) «V(1)*I(1) + V(2)*I(2) + V(3)*I(3) ; 
if ( t(k) >- 1/60 ) 

S » 0 ; 
for j * l:Th 

S - S + pl(k“j+l) ; 

end 

Pl(k) - S/Th ; 
elseif ( t(k) < 1/60 ) 

Pl(k) - 0 ; 

end 

XI - 1(1) - Pl(k)*V{l)/NormV ; 
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X2 « 1(2) - Pl(k)*V(2)/NormV ; 
X3 - 1(3) ~ Pl(k)*V(3)/NormV ; 
Current * [ XI ; X2 ; X3 ] ; 


A.3 Program for compensation of arc farnace with a current source 
realized through a 3phase 2IeveI inverter 


clear all 


global R1 R2 R3 R4 R5 Ro Rf Lf 
global Rs Is Rtl Ltl Rt2 Lt2 Ra Rb Rc 
global ja Ka jb Kb jc Kc 
global Fla Gla Fib Gib Flc Glc 

global Ala A2a Bla B2a Alb Bib A2b B2b Ale Blc A2c B2c 

global Xa Xb Xc 

global II Fc Iq h pi Th PI 

global Ic t Td Tc w Vs Vo Da Db Uc Vp Vdo Vp Vc Cf Ploss 
TC - 0.03 ; Tf - 1.5 ; Td - 2/(360*60) ; 

f - 60 ; w - 2*pi*£ ; Vo - 33000*sqrt (2/3) ; Vdc - Vo*1.2 ; Kp - 10000.0 ; 
h - 10 ; Th - round(l/(2*Td*f) ) ; Rf - 0.01 ; Lf - 0.015 ; Cf - 0.01 ; 

R1 - 3 ; R2 - 8 ; R3 - 10 ; R4 - 15 ; R5 - 20 ; 

Rs - 0.254 ; Ls » 0.00005 ; 

Rtl- 0.06 ; Ltl= 0.001 ; 

Rt2- 0.001 ; Lt2- 0.03 


)c -1 ; t()c) =0 ; Vc(l) = Vdc ; 

Vs(l,)c) - Vo*sin(w*t()c) ) ; 

Vs(2,k) = Vo*sin(w*t()c) - 2*pi/3) ; 

Vs(3,)c) - Vo*sin(w*t(k) + 2*pi/3) ; 
Xa(:,k)-[0;0]; 

Xb(:,k)-[0;0]; 

Xc(:,k) = [0;0]; 
Il(:,k)-[0;0;0]; 

Is ( : ,k) = [ 0 ; 0 ; 0 ] ; 

Iq(:,k) - [ 0 ; 0 ; 0 ] ; 

pl(l) ” 0 ; Pl(l) “ 0 ;ps(l) = 0 ; 

Fc(:,k) - [ 0 ; 0 ; 0 ] ; 

Ql(:,k) - [ 0 ; 0 ; 0 ] ; 

Qs(:,k)-[0;0;0]; 
ja ” 0 ; jb " 0 ; jc -■ 0 ; 
Ka-0;Kb-0;Kc-0; 

Ra » 0 ; Rb - 0 ; Rc - 0 ; 

Ua(:,k) - t V3(l,k) ; Fc(l,k) ) ; 

Ub(:,k) - [ V3(2,k) i Fc(2,k) ] ; 

Dc(:,k) - [ Vs(3,k) ; Fc(3,k) ] ; Coun = 0 ; 


while ( t(k) < Tf ) 


k - k + 1 ; 

t(k) = t(k-l) + Td ; 


solvejphaseA ; 
solve_phaseB ; 
5olve_phaseC ; 


11(1, k) 
11(2, k) 
11(3, k) 
ls(l,k) 
ls(2,k) 
ls(3,k) 


Xa(l,k) + Xa(2,k) ; 
Xb(l,k) + xb(2,k) ; 
Xc(l,k) + Xc(2,k) ; 
Xa(l,k) ; 

Xb{l,k) ; 

Xc(l,k) ; 
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Vp(l,k) - Ala*Xa{l, (k-1) ) + A2a*Xa (2, (k-1) ) + Bla*Vs (1, (k-1) ) + B2a*Fc(l, (k-1) ) ; 

vp(2,k) " Alb*Xb(l, (k-D) 4 A2b*Xb(2, (k-1) ) + Blb*V8 (2, (k-1) ) + B2b*Fc(2, (k-1 ) ) t 

Vp(3,k) - Alc*xc(l, (k-1) ) 4- A2c*Xc(2, (k-1) ) 4- Blc*Vs (3, (k-1) ) + B2c*Fc(3, (k-1) ) ; 

Vs(l,k) - Vo*sin(w*t (k) ) ; 

Vs(2,k) - Vo*sin(w*t (k) - 2*pi/3) / 

Vs(3,k) - Vo*sin(w*t (k) 4- 2*pi/3) ; 

updatecapV ; 

Ploas - Kp*(Vdc - Vc(k))/3 ; 

Coun • Coun 4- 1 ; 

Ql(:,k) - updateQ( Vp(:,k) , Il(:,k) ) ; 

Q 3 (:,k) - updateQ( Vs ( : , k) , Is(:,k) ) ; 

lq(:,k) - updatelq( Vp(:,k) , Il(:,k) ) ; 

ps(k) - Vs(l,k)*ls (l,k) 4- Vs(2,k)*ls(2,k) 4- Vs (3, k) *ls (3, k) ; 

if ( t(k) >“ 1/60 ) 

S-0; 
for j - l:Th 

S - S 4- ps(k-j4-l) ; 

end 

Ps(k) * S/Th ; 

end 

if ( t(k) > Tc ) 

if ( Xa(2,k) <- ( Iq(l,k) - h ) ) 

Fc(l,k) - Vc(k) - Vp(l,k) ; 
elseif ( Xa(2,k) >- ( Iq(l,k) + h ) ) 

Fc(l,k) - - Vc(k) - Vp(l,k) ; 
else Fc(l,k) - Fed, (k-D) ; 
end 

if ( Xb(2,k) <- ( Iq(2,k) - h ) ) 

Fc(2,k) - Vc(k) - Vp(2,k) ; 
elseif ( Xb(2,k) >- ( Iq(2,k) 4- h ) ) 

Fc(2,k) - - Vc(k) - Vp(2,k) ; 
else Fc(2,k) = Fc{2, (k-D) ; 
end 

if ( Xc(2,k) <- ( Iq(3,k) - h ) ) 

Fc(3,k) - Vc(k) - Vp(3,k) / 
elnalf ( Xc(2,k) >- ( Iq(3»k) 4- h ) ) 

Fc(3,k) - - Vc(k) - Vp(3,k) ; 
else Fc(3,k) - Fc(3, (k-D) ; 
end 

else Fc{:,k) « [ 0 ; 0 ; 0 ] ; 
end 

Ua(:,k) - [ Vs(l,k) ; Fc(l,k) ] ; 

Ob(:,k) - [ Vs(2,k) ; Fc(2,k) ] ; 

Uc(:,k) - [ Vs(3,k) ; Fc(3,k) ] ; 

end 



function updatecapV ; 


global k Cf Vc Td Xa Xb Xc 

A - - ( Xa(2,k) 4- xb(2,k) 4- xc(2,k) )/Cf ; 
B - 0 ; 

[F G] - c2d( A , B , Td ) ; 


vc(k) - vc(k-l) + A*Td ; 



> All the other functions being called in A.3 are the same as those for A.2 . 



